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ABSTRACT 

We have detected coherent oscillations ("dwarf nova oscillations") in Hubble Space 
Telescope spectra of the dwarf nova OY Car. The oscillations were seen towards the 
end of a superoutburst of OY Car. The oscillations are extraordinary compared to 
the many other examples in the literature for two reasons. First, their amplitude is 
large, with a peak-to-peak variation of 8 to 20% of the total flux over the range 1100 
to 2500A. However, most remarkably we find that there are two components present 
simultaneously. Both have periods close to 18 s (equivalent to 4800 cycles/day) but 
they are separated by 57.7 ± 0.5 cycles/day. The lower frequency component of the 
pair has a strong second harmonic while its companion, which has about twice its am- 
plitude, does not. The oscillation spectra appear hotter than the mean spectrum and 
approximately follow the continuum distribution of a black-body with a temperature 
in the range 30,000 to 50,000 K. 

We tentatively suggest that the weaker non-sinusoidal component could represent 
the rotation of the white dwarf, although we have been unable to recover any such 
signal in quiescent data. 
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1 INTRODUCTION 

Dwarf nova oscillations (DNOs) are one of the unsolved 
mysteries of cataclysmic variable stars. First discovered by 
Warner & Robinson (1972), DNOs appear during outbursts 
as moderately coherent oscillations (with Q — P/5P from 
10 4 to 10 6 ) with periods from around 7 to 40 seconds, (al- 
though recently the lower limit has been extended to 2.8 
seconds from observations of SS Cygni, van Teeseling 1997, 
Mauche 1997). The short periods clearly implicate the inner 
accretion disc and white dwarf. However while the rotation 
of the white dwarf is too coherent to match period changes 
seen in DNOs, it is hard to see how the accretion disc can 
produce anything as coherent as observed. The amplitudes 
(peak-to-peak) of DNOs in optical observations are typi- 
cally less than 0.5% and therefore they can often only be 
detected after period analysis. The amplitudes increase to- 
wards shorter wavelengths and are tens of percents in X-ray 
light curves (Cordova et al. 1984). The periods of DNOs 
correlate with system brightness, becoming shorter as the 
system becomes brighter and vice versa. 

DNOs seen so far have exhibited a single, sinusoidal 
signal (although quasi-periodic oscillations are occasion- 
ally seen as well). In this paper we present data that vio- 
lates both generalisations. We will show that DNOs present 
in Hubble Space Telescope observations of the short-period 



(P = 91 min) eclipsing dwarf-nova OY Car have two compo- 
nents, one of which has a strong second harmonic. We begin 
by describing the observations. 



2 OBSERVATIONS 

On 24 April 1992 (pre-COSTAR) we took 419 spectra cov- 
ering the range 1150 to 2510A with the G160L grating and 
the Faint Object Spectrograph (FOS) on the Hubble Space 
Telescope. Each exposure was 4.74 s long, and the time from 
the start of one exposure to the start of the next was 5.60 s. 

With the G160L grating, FOS also captures the zeroth 
order undispersed light which is weighted to longer wave- 
lengths. The pass-band of the zeroth order light has been 
calibrated by Eracleous et al. (1994) who found that it has 
a full width at half maximum of 1900A centred on 3400A 
and who determined a scale factor of 820 counts/sec/mjy 
for our 4.3" aperture. 

OY Car at this time was right at the end of a super- 
outburst which started some 17 days earlier. All amateur 
measurements taken after our run were upper limits, and 
so we believe that our data were taken within a day of the 
return to quiescence. Eight other HST observations of the 
outburst were made prior to those we describe here, the 
closest being observed two days earlier. No oscillations were 
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Figure 2. The light curve of OY Car in the first order dispersed light is plotted in the left panel whilst the right panel shows the zeroth 
order light curve. Both light curves show periodic oscillations. The oscillation is relatively stronger in the shorter wavelength (left-hand) 
panel. 
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Figure 1. The mean spectrum of OY Car observed at the end 
of a superoutburst with the G160L grating and HST/FOS. The 
lower line shows the la uncertainties on the spectrum. 



found in any of these. Similarly, we find no oscillations in 
any of 7 succeeding observations during quiescence, the first 
of which took place 10 days after the observations reported 
here. 



3 RESULTS 

The mean spectrum during the observations is presented in 
Fig. [lj. The spectrum is blue (note that it is plotted in terms 
of /„) with modest CIV emission. The continuum is very 
complex with many features from what has previously been 
recognised as the "iron curtain" of material that partially 
absorbs light from the inner disc and white dwarf (Home et 
al. 1994). The broad dips around 1600 to 1800A and beyond 
2200A are highly characteristic in this regard. 

It is the light curves (Fig. |^) that prove to be unusual. 
As well as very deep eclipses characteristic of a small light 
source concentrated at the centre of the disc, they exhibit 
what appears at first to be a curious noise pattern. This is in 
fact the dwarf nova oscillation, and is a rare case in which it 
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Figure 3. The Lomb-Scargle periodograms of the two orders are 
plotted beyond the Nyquist frequency (dashed line). The peak 
at 5900 cycles/day is an alias of the second harmonic (seen near 
9500 cycles/day) of the component near 4800 cycles/day. 



can be seen directly in the raw data. The mean flux level of 
around 9 mjy compares to a level of about 1 mjy observed 
during quiescence. 

The periodogram (Lomb 1976; Scargle 1982) of each 
light curve is shown in Fig. |j| These were computed after 
masking out data affected by the eclipse and then fitting 
and subtracting a 7th order polynomial from the data to 
remove long timescale variations. At high frequencies, the 
periodograms are dominated by two peaks close to 4800 and 
5900 cycles/day respectively. We will show below that the 
lower frequency "peak" is composed of two closely spaced 
signals. After some time we realised that the higher fre- 
quency peak is the alias of a peak near 9500 cycles/day as we 
show in Fig. ^ by extending it beyond the Nyquist frequency 
(7719 cycles/day). We make this identification because, as 
we will show, the 9500 cycles/day peak is then the second 
harmonic of the weaker and lower frequency component of 
the 4800 peak. The aliasing is a consequence of the length of 
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Table 1. Frequencies of components 

Bayesian models 



4600 4B00 5000 9400 9600 9800 

Frequency (cycles/day) Frequency (cycles/day) 

Figure 4. The figure shows regions around the two main peaks 
of the periodogram of Fig. In the left-hand panels the peak 
marked with an arrow is stronger in the data than in the window 
function (generated by replacing the data with a pure sinusoid). 
On the right, the arrows indicate the expected position of the 
second harmonic of the peak marked by an arrow in the left-hand 
panels. There is no sign of a second harmonic of the stronger peak 
of the left-hand panel. 



our exposures, which have the additional effect of reducing 
the amplitude of the signal. If a signal of period P is sampled 
with exposures of length At, then the observed amplitude 
is sm(8)/8 times the true amplitude where 8 — nAt/P. In 
our case At — 4.74 s, whereas 9500 cycles/day corresponds 
to P = 9.09 s, leading to a reduction factor of 0.61. In the 
periodogram, this is squared and means that the peak at 
6000 cycles/day is only 37% of its true height. For the same 
reason, the third harmonic is reduced to only 7% of its true 
height. Unfortunately the aliasing means that the third har- 
monic is expected to be at 1160 cycles/day in a region of 
high background from low frequency power, and it cannot 
be detected. 

We said above that the main peak has two components 
the weaker of which has a second harmonic. This is shown 
in Fig. ^ where we isolate the two peaks at 4800 cycles/day 
and the second harmonic at 9500 cycles/day. We also dis- 
play "window" functions computed by replacing the data 
by pure sinusoids, including the mask around eclipse. In the 
left-hand series of panels, the data show a peak immediately 
to the left of the main peak which is not seen as strongly 
in the window function. The strongest peak is at 4816.32 
cycles/day while its companion is some 53 cycles/day lower 
at 4763.12. Unfortunately because of the short length of our 
observation, there is a weak peak at the same place in the 
window function and on this basis alone we would be hesi- 
tant in claiming the presence of two signals. However in the 
right-hand panels, we see the second harmonic of the 4763 
component while there is no detectable second harmonic of 
its stronger companion. This shows that the component in- 
dicated by an arrow in the left-hand panels of Fig. ^ is not 
merely an artefact of sampling. 

The Lomb-Scargle periodogram is not really appropri- 
ate in cases where more than one frequency is present, al- 
though if they are far enough apart we would not expect 
significant problems in determining the frequencies of the 
peaks. In this case however, two frequencies are close, and 
there must be doubt over the accuracy of the numbers 
quoted above. On the other hand, we need to know them 



Lomb- 
Scargle^ 
4816.32 
4763.12 
9518.41 



1 cpt 

4816.26 ± 0.27 
4763.19 ± 0.42 
9518.36 ± 0.65 



3 cpt 
4816.59 ±0.34 
4760.13 ± 1.21 
9517.94 ±0.68 



2 cpt (/ 3 =2/ 2 ) 

4816.72 ±0.31 
4759.05 ± 0.33 
9518.10 ±0.66 



The frequencies are all measured in units of cycles/day 



accurately in order to verify quantitatively our identification 
of the second harmonic. That this is a significant problem is 
seen once uncertainties are placed on the frequencies. Using 
Bayesian probability theory (as outlined in the appendix) 
it is straightforward to extend the periodogram to account 
for any number of periodic components. For a single compo- 
nent, this reduces to a form very close to the Lomb-Scargle 
version, and shows that the periodogram is closely related 
to the natural log of the posterior probability distribution 
of the period. The curvature then gives the uncertainty, and 
applying this we obtain the frequencies of the three main 
peaks listed in the second column of Table which can be 
compared to the values obtained from the Lomb-Scargle pe- 
riodogram listed in the first column. As expected the values 
in the first two columns of table ji] are almost identical, but 
now with uncertainties we can compare the frequencies of 
the peaks at / 2 = 4763 and /3 = 9518 cycles/day that we 
claim are harmonically related to each other. We find the 
difference 2/2 — fa — 8.02 ± 1.06 cycles/day, and it would 
seem that there is a problem. This is the reason why in 
Fig. ^ the predicted second harmonic indicated by the arrow 
in the right-hand panels is not a perfect match to the ob- 
served peak. However, once three components are included 
correctly (see the appendix), we obtain the values in the 
third column. The frequencies do indeed change, and now 
the difference between the harmonically related components 
becomes 2/2 — /a = 2.32 ± 1.39, in acceptable agreement 
within the uncertainties. Given the harmonic relation, we 
can re-fit the frequencies on the basis that /j = 2/2 exactly, 
so that there are only 2 independent frequencies, and we 
arrive at the fourth column which contains the final values 
that we will use from now on. It can be seen that it is the 
harmonic at 9518 cycles/day which dominates the determi- 
nation of the frequency, and as we will show, it has a higher 
amplitude than the fundamental. 

The referee has pointed out that a single component 
with a phase and/or period change might appear as more 
than one period in periodograms, and that just such changes 
have been observed in other examples of DNOs. This is cer- 
tainly true. We carried out some tests with a single sig- 
nal of constant period suffering a sharp phase shift about 
half way through the observations. As the phase shift in- 
creases, one of the two side peaks grows in strength while the 
main peak shifts in frequency and weakens. When the phase 
shift reaches about 120° the periodogram appears not unlike 
Fig. ^j. In this case the "true" frequency is placed somewhere 
between the two peaks, but is closest to the strongest peak. 
The coincidence of the weaker peak (the /2 component) with 
one of the two side-lobes of the /1 component is suspicious 
in this regard. We reject this possibility however because of 
the /3 component which agrees to within 1.7c of 2/2 and 
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Figure 5. The phases of the 3 components of the oscillations are 
plotted versus wavelength. The right-most point is the phase of 
the zeroth order light curve. There is no evidence for any depen- 
dence of phase upon wavelength. The phases are measured with 
respect to the zero point HJD = 2448736.835. 



yet is 120<t away from 2/i. This means that fz is the second 
harmonic of f% , and thus component fa is genuinely indepen- 
dent of fi. We further note that while we know of no other 
multiply periodic DNOs, multiple periodicity has been seen 
in the dwarf nova WZ Sge in quiescence (Patterson et al. 
1998). 



3.1 Spectra of the oscillations 

Armed with accurate frequencies for the three components 
(two of them harmonically related), we can fit sinusoids at 
many different wavelengths to obtain the spectra of the oscil- 
lations. To do so we fit all three components simultaneously 
to the spectra, using the form a cos 2irft + b sin 2n ft for each 
component. This can be re-expressed as Acos2Tv(ft + (j>) 
where the semi-amplitude A — \/a 2 + b 2 and the phase <j> 
is given by t&n2n<f> — —b/a. The estimate of amplitude A 
suffers from a noise bias (it can only be positive). Although 
this can be corrected approximately, it is preferable to see 
first if there is any evidence for a variation of phase with 
wavelength. If not, then <f) can be fixed and we can derive an 
unbiassed estimate of the semi-amplitude. The phase versus 
wavelength is plotted in Fig. ^ and shows that the phase is 
indeed independent of wavelength in all three components. 
The weighted mean phases of both orders was then com- 
puted in order to fit the amplitudes. We obtain the values 
110.6° ± 0.8, -15.4° ± 1.9 and -57.7° ± 1.9 for the 4816, 
4759 and 9518 components respectively. The latter two val- 
ues show that the 9518 harmonic peaks 13.9° ± 2.1 after the 
4759 fundamental, as measured in terms of the fundamen- 
tal's phase. 

Holding the phases fixed we then fitted the amplitudes 
which are plotted Fig. [| The oscillation spectra show no sign 
of the CIV emission seen in the mean spectrum (at 1550A). 
The CIV is thought to form by scattering of continuum pho- 
tons in an outflow above and below the disc. It therefore 
gives us an angle-integrated view of the inner disc bright- 
ness. This suggests that the oscillations are not produced 
by the entire luminosity of the system varying, but rather 
by the varying visibility of some hot source. This counts 
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Figure 6. The semi-amplitudes of the 3 components of the oscil- 
lations are plotted versus wavelength along with the mean spec- 
trum binned to the same scale. The right-most point is the zeroth 
order point scaled according to the calibration of Eracleous et al. 
(1994) and shifted in wavelength to avoid compressing the first 
order spectrum; while the statistical uncertainties on these points 
are low there is a 50% uncertainty in their calibration which would 
move them all up or down by the same factor. The dashed lines 
show black-body spectra for T = 30,000 K, while the dotted lines 
are Raleigh-Jeans spectra. The values plotted are not corrected 
for the amplitude reduction caused by the finite length exposures. 



against radial oscillation models such as those of Molteni et 
al. (1996). 

The oscillation spectra are very much hotter than the 
average spectrum. This can be seen over the range of the 
first order dispersed light but more clearly still in the ratio 
of first to zeroth order light. This agrees with earlier work 
at optical wavelengths (e.g. Middleditch & Cordova 1982; 
Stiening et al. 1984); our data extend this result into the 
space ultraviolet. Partly as a result, the peak-to-peak am- 
plitude in the first order light is very high, reaching almost 
20% at the shortest wavelengths of the 4816 component. 
In the zeroth order light, this has fallen to 5%. Schoembs 
(1986) observed a maximum amplitude of 2.2%, but given 
that his observations were approximately B band (our ze- 
roth order has a central wavelength of 3400A) and the very 
blue spectra, our values are probably comparable. 

The spectra of DNO are so blue that ground-based data 
can only put a lower limit upon their temperature. With our 
coverage of the ultraviolet we might hope to fare better. In- 
deed, the Raleigh- Jeans spectra plotted in Fig. ^ are hotter 
than any of the oscillation spectra, particularly in the small 
amount of zeroth order flux they produce. Allowing for the 
50% uncertainty Eracleous et al. (1994) quote on their cal- 
ibration of the zeroth order flux, we deduce an upper limit 
of about 50,000 K if the oscillation spectra are black-bodies. 
The lower limit is independent of the zeroth order calibra- 
tion since it is fixed by the short wavelength ultraviolet flux. 
We place a lower limit of about 30,000 K (see Fig. ^) from 
such a comparison. 

The ratios of the oscillation amplitude in the first order 
light divided by its amplitude in zeroth order (both mea- 
sured in counts/sec) are 1.37T0.04, 1.10T0.08 and 1.53T0.11 
for the 4816, 4759 and 9518 components respectively; the ra- 
tio of the mean fluxes in each order is 0.68. Thus the second 
harmonic 9518 component is hotter than its fundamental 
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Figure 7. The shape of the 4759/9518 obtained from sinusoid 
fits and corrected for the amplitude reduction caused by the finite 
exposure lengths. The solid and dashed lines represent the first 
and zeroth orders respectively. The pulse phase is measured in 
terms of the 4759 fundamental with respect to the zero point 
HJD = 2448736.835. 



counterpart at 4759 cycles/day. Moreover, the amplitude of 
the second harmonic is larger than the fundamental. We find 
the ratios of the 9518 to 4759 amplitudes (corrected for the 
finite exposure factor discussed earlier) to be 1.58 ± 0.09 in 
the first order and 1.14±0.09 in the zeroth order. Along with 
the phases mentioned earlier, these lead to the pulse shapes 
shown in Fig. |^. This plot shows a significant "interpulse" 
which will be of significance in the discussion section. The 
true pulse shape may be more complex since our exposures 
were too long to detect any but the first two harmonics. 



3.2 The eclipse of the oscillations 

As OY Car is eclipsing, we can learn something about the 
distribution of the oscillating light by measuring the ampli- 
tude and phase of the oscillations as a function of time. The 
standard method of doing so is to hold the period fixed and 
then fit sinusoids to short sections of the data. How short a 
section can be used depends upon the strength of the oscil- 
lations. Our oscillations are especially strong, but we have 
the extra complication of multiple periods. While it is al- 
most as easy to fit to several periods as it is to fit to one, 
the solutions become ill-conditioned if data segments are not 
long enough to separate the periods. This is most obviously 
a problem with the 4816 and 4759 components which can 
only barely be separated using all of our data. Therefore 
for the purposes of this section we simply ignore the weaker 
component (4759) while realising that it will cause a certain 
amount of interference or "beating" . Unfortunately spectral 
leakage from the 4816 component can also affect the 9518 
component, which since it appears at 5900 cycles/day be- 
cause of aliasing is much closer than it seems. We applied a 
linear taper to the data segments in an effort to reduce the 
leakage, but were forced to use segments of 30 points each 
for the 9518 component. For the 4816/4759 combination we 
settled on 16 points per segment, and in each case we used a 
50% overlap from one segment to the next. Prior to splitting 
the data up we also applied a high-pass filter to remove low 
frequencies in order to prevent leakage from them affecting 



the results. The results of the fits are shown in Figs, g and 
for the first and zeroth order light respectively. 

In interpreting the plots it should be remembered that 
the amplitude measurements are subject to a positive bias 
of the same order as the uncertainties plotted and so the 
amplitudes during eclipse are consistent with no detection 
even though it appears that with several marginal detections 
there is some evidence for oscillations during eclipse. Both 
the mean flux and oscillations are deeply eclipsed, but from 
the zeroth order results on the 4816 component in Fig. |we 
can say that the oscillations are even more deeply eclipsed 
than is the mean light. The other panels give no clear answer, 
partly because the eclipse is very deep even in the mean 
light, and partly because of the noise bias discussed earlier. 

Since we know that the period of the oscillations imply 
an origin in the inner disc or white dwarf, it may not seem 
that we have learned much. However, in other systems, such 
as DQ Her and UX UMa, the oscillating light is spread out, 
indicating that we see light reprocessed from the disc. Phase 
shifts observed during eclipse show where the reprocessing 
occurs (Petterson 1980). It seems then that such effects play 
at best a minor role in ou r da ta; we will look at reprocessing 
in more detail in section 4.1. 



4 DISCUSSION 

Our HST data are qualitatively different from previous ob- 
servations of dwarf nova oscillations. The essential features 
of our data are (a) the simultaneous presence of two oscil- 
lations each with a period close to 18 s but separated by 
57.7 ±0.5 cycles/day, and (b) a strong second harmonic of 
the weaker and lower frequency component which peaks just 
13° after the fundamental, and no second harmonic of the 
stronger component. 

To our knowledge, there are no previous cases of either 
multiple or non-sinusoidal oscillations. As has long been re- 
alised, the short periods of DNO imply an origin in or near 
the white dwarf. The direct interpretation of our data is that 
this region is capable of producing two oscillations of differ- 
ent character at the same time. Although it is not clear why 
it has not been seen before, it is what we believe is happen- 
ing. However, before proceeding upon this basis, we will first 
consider an attractive alternative which turns out not to be 
satisfactory. 

4.1 Reprocessing models 

We already know of one class of cataclysmic variable star in 
which a signal at one period can give rise to others. These 
are the DQ Her stars in which a coherent signal from the 
rotation of a magnetic white dwarf can be detected. These 
stars also show sidebands to the rotation frequency displaced 
by the orbital frequency. Such sidebands can result from 
modulation of the spin frequency u) on the orbital frequency 
Q which leads to u)±Q sidebands. A single sideband at cj— 
can result from reprocessing of the spin pulses off structures 
fixed in the rotating frame or from direct interaction of the 
gas stream with the magnetosphere of the white dwarf. It is 
clear then that such effects are possible, but could they be 
occurring in OY Car? 

The DQ Her analogy suggests that the 4759 cycles/day 



© 1998 RAS, MNRAS 000, 000-000 



6 T. R. Marsh, Keith Home 



4816 



9518 



o 

m 



o 
o 



o 
m 



o 
oo 

CO 

N 
O 

o 



o 

CO 



o 
o 

CO 



r 

- 


i 1 

i •••• 

k 

4]h — 1 — 1 — | — 1 — h- 


■ 

i | i i i i | i i i \'"\" 


- 1 T fi * 

+ 

1 


i vv 

+j-H — i — i — | — i — h- 


1 1 1— 1 1 1 1 1 1 1 jf"|«» 

-iW V ++ + + : 

~ * ' T II 


i : 

+|h — i — i — | — i — i— 


— i — i — |— i — i — i — i — — i — i — i il^ !■ 
••••••••• • 1 

... .... ......... 

••••••••• • ' 

, i , * , i * in j 


• l _ 

• 1 _ 



-0.2 



-0.1 



0.1 



o 

m 



o 
o 



o 
in 



o 

CO 



o 
o 
co 

o 
o 

CO 



f-^ — I — I — I — I — I — I— I * I — I — I — I — I — 



\ + 1 

I I 1 1 I . I . i I . 1 . I 



-0.2 



-0.1 



0.1 



Orbital phase - 75146 



Figure 8. From the top are plotted the mean flux, oscillation amplitude, percentage oscillation amplitude and oscillation phase (degrees) 
evaluated over data segments of 16 points for the 4816 component and 30 for the 9518 component for the first order light. 



component might be a sideband of the 4816 component. We 
expect it to be this way round because the 4816 compo- 
nent is the stronger and because if the light source and re- 
processing structure orbit in the same direction, the repro- 
cessed frequency is lower than the main signal. In addition 
the 4816 component is sinusoidal and thus can be identified 
with the "standard" DNO of other observations. The repro- 
cessing structure must advance in the inertial frame at the 
57.7 ± 0.5 cycles/day splitting we measured earlier. This is 
3.6 ± 0.03 times the orbital frequency in OY Car. During 
superoutburst it is thought that the accretion disc reaches 
a radius where it orbits 3 times for each orbit of the binary 
(Whitehurst 1988). The 3:1 resonance then drives the disc 
into an eccentric shape and subsequent precession can then 
explain the observation of flares or "superhumps" on periods 
longer than the binary period. From Kepler's third law, a 
ratio of 3.6 corresponds to a radius of about 89% of the 3:1 
radius. Since our observations were placed right at the end 
of outburst when the disc should have shrunk somewhat, 
this seems plausible. There are many other constraints that 
such a model must satisfy however. 

First of all there must be an asymmetry in the outer 
disc so that we can see it lit up every time the beacon from 
the inner disc sweeps by it. One cannot expect any such 
asymmetry to last very long since shear in the disc always 
tends to smear out azimuthal structure. While we cannot 
suggest a mechanism for generating such an asymmetry, its 
short-lived nature could be viewed as a plus point since then 
such extra components would be rarely seen. In order to be 



consistent with the eclipse results, we also have to suppose 
that the asymmetry is small enough to be totally eclipsed. 

Even allowing the above constraints, there are higher 
hurdles for reprocessing models to face. The blue colours of 
the oscillation spectra imply a high effective temperature, 
particularly so since our spectra extend into the space ul- 
traviolet. A comparison of the zeroth to first order flux ra- 
tio with black-body models of reprocessing suggest an effec- 
tive temperature of at least 20,000 K. This is lower than the 
30,000 K limit discussed earlier because oscillation spectra 
can appear more blue than the source of light that produces 
them since they depend upon the derivative of the spec- 
trum with respect to temperature. The luminosity L needed 
to produce reprocessed radiation with effective temperature 
T e ff at radius R in the disc is given by 

L = A^oT^- dXnR 



d(H/R) ' 

where H is the scale height in the disc. Taking H/R oc R 13 , 
then the derivative term becomes l/(/3H/R). Typically /3 ~ 
1/8 (Shakura & Sunyacv 1973) and we take H/R ~ 1/20 
and so the derivative term is of order 160. 

Taking Wood & Home's (1990) parameters for OY Car, 
the 3.6:1 radius corresponds to R — 0.41a or 1.9 x 10 8 m. 
We then find that to reach our minimum of 20,000 K the 
central peak luminosity must be of order L — 1700 Lq. The 
accretion luminosity of the inner disc and boundary layer, 
GM M / R, set an upper limit to the mean luminosity of the 
oscillating source (which is half the peak luminosity assum- 
ing 100% modulation). We find M > 3 x 10~ 7 M /yr, which 
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Figure 9. From the top are plotted the mean flux, oscillation amplitude, percentage oscillation amplitude and oscillation phase (degrees) 
evaluated over data segments of 16 points for the 4816 component and 30 for the 9518 component for the zeroth order light. 



is higher than the maximum rate expected for dwarf novae, 
and a strong point against reprocessing. 

The second harmonic is perhaps even more difficult to 
understand on the reprocessing model. It requires a non- 
linear repose to the incident flux and can be generated if as 
the flux incident on the disc rises, the reprocessing efficiency 
rises as well, thus sharpening the peaks. However even if 
this process is taken to its extreme to produce a series of 
sharp spikes or delta-functions, we can only obtain a second 
harmonic equal to the fundamental. In our case the second 
harmonic is stronger than the fundamental by a factor of 
1.58±0.09 in the first order light and 1.14±0.09 in the zeroth 
order as illustrated in Fig. []. In order to obtain the pulse 
shapes seen one is forced to suppose a complex response 
in which the reprocessing efficiency at first decreases and 
then increases with rising incident flux. We find this very 
implausible. 

Our conclusion is that a single oscillation plus repro- 
cessing cannot explain our observations: there really are two 
oscillations present at the same time. 

4.2 Two oscillations 

Accepting that both oscillations are generated on or near 
the white dwarf, how do they originate? The trapped gl- 
and r-mode model of Papaloizou & Pringle (1978) can cer- 
tainly produce more than one oscillation frequency. Indeed, 
it is hard to see how this model could only produce one fre- 
quency, an argument that has been used in the past against 
it. Now we have seen more than one component. However, in 



our view this does not support Papaloizou & Pringle's model 
because each component has a very different character, and 
if two oscillations are possible, why not more? 

Another model of DNOs frequently discussed is one 
in which the DNO frequency represents the keplerian or- 
bital frequency of some part of the inner disc (Bath 1973; 
Mauche 1996). On this model for instance, the 4816 com- 
ponent might represent the keplerian frequency, and, if so, 
it corresponds to an orbital radius of 1.4 times the radius 
of the white dwarf on the physical parameters of Wood & 
Home (1990). Explaining the relatively high coherence of 
DNOs has always seemed difficult for keplerian models as 
why should just one specific orbital frequency dominate? 
Our data perhaps provide some support however because 
there is an appealing explanation for the 4759/9518 com- 
bination, which might represent either the white dwarf ro- 
tation frequency or the beat frequency between the white 
dwarf and keplerian frequencies. In the latter case the white 
dwarf would rotate at the frequency difference of 57.7 ± 0.5 
cycles/day. An attractive feature of this model is that it is 
entirely natural that the two components should have a dif- 
ferent nature and second harmonics are easily produced if 
the white dwarf has two accreting poles for instance. Given 
a choice between 4759 or 57.7 cycles/day for the rotation 
frequency of the white dwarf, we would favour 4759 since 
we already know of short rotation period examples such as 
AE Aqr and DQ Her, amongst weakly magnetic accretors, 
a category which OY Car surely falls into. On the contrary, 
it is hard to imagine how the white dwarf in OY Car could 
have the 25 minute period required to match the 57.7 fre- 
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quency splitting as this would be longer than in many of the 
longer period DQ Her stars. On the basis of this idea, we 
searched for coherent periodicities in 7 epochs of HST taken 
during quiescence following the outburst and in Wood et 
al.'s (1989) ground-based photometry. We failed to find any 
coherent signals, with detection levels of 0.1 and 0.3% re- 
spectively (these limits only apply above 1000 cycles/day, 
because orbital variations and nickering limit the sensitivity 
at low frequencies). In this respect at least the 4759/9518 
component behaves like other DNOs which have never been 
seen in quiescence. 

It is not clear how much the failure to detect the 
4759/9518 component in quiescence counts against its be- 
ing the rotation frequency of the white dwarf. We speculate 
that perhaps it only becomes evident when the inner disc or- 
bital period (represented by the 4816 component) is close to 
co-rotation with the white dwarf (4759 component). Such a 
situation may favour the largest difference between accretion 
rate onto the poles compared with other parts of the white 
dwarf making the white dwarf spin period detectable. We 
note that the shortest DNO period from OY Car observed 
by Schoembs (1986) of 19.44 seconds is over 300 cycles/day 
away from 4759 and may have already been too different to 
generate the 4759/9518 signal. 

A third, hybrid model has been put forward by Warner 
(1995), based upon an earlier model of Paczynski (1978). In 
this model the DNO frequency is set by the outer layers of 
a weakly magnetic white dwarf decoupling from the white 
dwarf itself. The variation in period is put down to the spin 
up/down of the layer with accretion rate. This model suffers 
from the same problem as the keplerian models in that it is 
hard to see how a coherent oscillation can result. Now with 
our data, there is an additional difficulty as Warner's model 
does not appear to have any natural means for producing 
more than one periodicity. 

In summary while none of the current models provide 
a convincing explanation for DNOs, our observation of two 
components seems more naturally explained by keplerian 
orbit models than any other. 



5 CONCLUSIONS 

We have found dwarf nova oscillations at the end of a super- 
outburst observed with HST. The oscillations have a large 
amplitude, reaching 20% peak to peak at the shortest ul- 
traviolet wavelengths. Most remarkably there are two os- 
cillations present with similar periods, and the weaker and 
longer period of the two has a strong second harmonic. 

We discuss several current models and while none are 
persuasive, models in which the DNO frequency represents 
an orbital frequency of the inner disc give the most natural 
explanation of the additional component as it can represent 
the rotation of the white dwarf. However, we have been un- 
able to find any sign of such a component during quiescence, 
and so, if our suggestion is true, some means of suppressing 
the signal is required. 
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APPENDIX A: BAYESIAN PERIOD 
MEASUREMENT 

The detection and measurement of periodic components can 
be handled using Bayesian probability theory (e.g. Gregory 
& Loredo 1992; Bretthorst, 1988). We consider the proba- 
bility that our data is consistent with zero (model Z) versus 
the probability that it can be represented by a sum of N 
sinusoids (model S). In practice, these models imply that 
data must be "de-trended" by fitting polynomials or splines 
to remove low frequency variability prior to analysis. 

We wish to compute the probability ratio of these mod- 
els given the data, P(S\D) / P{Z\D). Let the flux at time U 
be yi, with uncertainty cr», for i — 1 to M. Our sinusoid can 
be represented in linear form as 
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Vi = 



2N 
3=1 



where the Xij functions are alternately the value of cosine or 
sine evaluated for frequency j at i*. Thus Xi t i = cos2nfiti, 
%i,2 = sin27r/iti, Xi^ = cos2nf2ti, etc. Written in this fash- 
ion rather than with phases makes later integrations simpler. 
B ayes' theorem gives 

P(S\D) _ P(S) P(D\S) 
P(Z\D) ~ P(Z) P(D\Z)' 

The first term on the right expresses one's prior belief in 
each model. It is important when one is discussing whether 
a component is really present, but since it remains constant, 
it does affect the periods measured. Information about the 
periods is contained in the second term, which is the ratio of 
the probability of our obtaining the data given the models. 

Let us first consider P(D\Z). We assume that the data 
are gaussian random variables then 



P{D\Z) = 



1 



II 



exp 




The probability of the data given the sinusoid model is more 
complicated because we first need to compute it for a spe- 
cific vector of parameters, i.e. P(D\S, a). This is straight- 
forwardly 



P(D\S,a) = 



1 



exp ■ 



E 



2a? 



In order to obtain the probability we want, P(D\S), we 
must integrate over the parameters a: 



P(D\S) = 



P(D\S,a)P(a) da. 



Once again P{a) is a prior probability, in this case of the 
semi- amplitudes. We take the amplitudes to have indepen- 
dent, uniform distributions over the range —R/2 to +R/2 
so 



Ei W i X i,2 



and the IN by IN symmetric matrix X is given by 



X = 



/ Ei w i x i,l x i,l 

E 4 mxi,2%i,i 



J2i W i X iA X i,2N \ 
Ei w i x i,2 x i,2N 

Ei 1Vi x i,2NXi,2N ) 



where Wi = l/of and the sums extend from i = 1 to M, the 
number of data points. 

Finally dividing by P(D\Z) and taking the natural log- 
arithm we obtain 



. P(D\S) nAri 
In ' { = 2N In 
P(D\Z) 



2jt 



R 



iln(detX) + ib*X _1 6. (Al) 



Together the first two factors represent the ratio of the vol- 
ume of parameter space allowed by the data compared to 
the volume (R 2N ) allowed by the prior model of the semi- 
amplitudes. This is a factor penalising the sinusoid model for 
its adjustable parameters. The last factor is the important 
one. For example, if we take a single sinusoid model with 
uniform uncertainty values, one can show that the value of 
b'jf - 6/2 is identical to the expression given by Scargle 
(1988) for a periodogram of unequally spaced data. 

The best frequencies are those which maximise the log 
probability ratio of Eq. A.1 (which is why R has no effect 



upon their values). The maximum can be located using stan- 
dard techniques. If at the same time the second derivative 
matrix is evaluated, then uncertainties can be computed. 
This is the procedure we applied to derive the frequencies 
listed in table of this paper. 



P(a) 



R 2N- 



One normally has little idea of what R would be, so it should 
be set conservatively. For example, setting R equal to the 
brightest OY Car has ever been seen would be a reasonable 
course to take. Its value has some effect on the significance 
of a peak but not on the value of the period deduced. 

With these assumptions we now need to integrate the 
exponential over a 2iV dimensional cube of side 7?. To make 
this simpler we assume that this cube contains all the region 
over which the exponential is significant. After some algebra 
we obtain 



P(D\S) 



exp 



R 




where the vector b is given by 
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